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The  dynamics  of  a  large  number  of  liquids  and  polymers  exhibit  scaling  properties  characteristic  of  a  simple 
repulsive  inverse  power-law  potential,  most  notably  the  superpositioning  of  relaxation  data  as  a  function  of  the 
variable  TVY ,  where  T  is  temperature,  V  the  specific  volume,  and  y  a  material  constant.  A  related  scaling  law 
TmVy,  with  the  same  exponent  T  =  y,  links  the  melting  temperature  Tm  and  volume  Vm  of  the  model  IPL 
liquid;  liquid  dynamics  is  then  invariant  at  the  melting  point.  Motivated  by  a  similar  invariance  of  dynamics 
experimentally  observed  at  transitions  of  liquid  crystals,  we  determine  dynamic  and  melting-point  scaling 
exponents  y  and  T  for  a  large  number  of  nonassociating  liquids.  Rigid,  spherical  molecules  containing  no 
polar  bonds  have  F  =  y;  consequently,  the  reduced  relaxation  time,  viscosity,  and  diffusion  coefficient  are  each 
constant  along  the  melting  line.  For  other  liquids  y  >  F  always;  that  is,  the  dynamics  is  more  sensitive  to  volume 
than  is  the  melting  point,  and  for  these  liquids  the  dynamics  at  the  melting  point  slows  down  with  increasing  Tm 
(that  is,  increasing  pressure). 

DOI:  10.1 103/PhysRevE.83.031504  PACS  number(s):  64.70.pm,  64.70.D- 


I.  INTRODUCTION 

When  cooled  or  compressed,  liquids  undergo  a  series  of 
transitions  associated  with  qualitative  changes  in  dynamical 
properties  [  1  ] .  At  a  relatively  high  temperature  or  low  pressure, 
molecular  cooperativity  of  the  motions  arises,  evidenced  by 
the  onset  of  non- Arrhenius  and  non-Debye  behaviors.  At 
lower  temperature  (or  higher  pressure)  a  dynamic  crossover 
is  observed,  associated  with  a  characteristic  change  in  the 
temperature  dependence  of  the  relaxation  time,  along  with 
deviations  from  the  Stokes-Einstein  and  Einstein-Debye  re¬ 
lations.  Eventually,  at  a  sufficiently  low  pressure-dependent 
temperature,  translational  and  rotational  motions  cease  as  the 
material  becomes  a  glass.  Studies  of  liquid  dynamics  as  a 
function  of  temperature  and  pressure  have  shown  that  for  a 
given  liquid  the  relaxation  time  is  constant  at  each  of  these 
three  dynamical  transitions:  the  onset  of  cooperativity,  the 
dynamic  crossover,  and  the  glass  transition.  That  is,  although 
the  transition  temperatures  (TA,TC,  and  Tg)  depend  on  pressure, 
the  corresponding  relaxation  times  at  the  transition  are,  to 
a  very  good  approximation,  independent  of  thermodynamic 
conditions  [1]. 

For  many  liquids,  these  transitions  cannot  all  be  observed, 
usually  due  to  crystallization,  a  genuine  thermodynamic  transi¬ 
tion.  For  liquid  crystals,  which  are  rigid  anisotropic  molecules, 
the  crystalline  solid  is  reached  via  a  series  of  transitions 
through  mesomorphic  (nematic,  smectic,  etc.)  states,  which 
exhibit  orientational  order  but  limited  or  nonexistent  long- 
range  translational  order.  Surprisingly,  although  the  transition 
temperatures  per  se  are  pressure-dependent,  the  relaxation 
time  for  longitudinal  reorientation  (around  the  short  axis  of 
the  molecule)  of  liquid  crystals  is  constant  at  the  transitions 
from  mesomorphic  states  to  the  isotropic  liquid  and  the 
crystalline  solid,  as  well  as  between  different  mesomorphic 
states  [2,3] .  This  constancy  of  r  has  thus  far  only  been  observed 
for  mesogenic  liquids.  Whether  the  phenomenon  extends  to 
simple,  isotropic  liquids  is  unknown. 

It  is  observed  experimentally  that  for  a  wide  variety  of 
substances,  including  van  der  Waals  liquids,  polymers,  and 


even  ionic  liquids  and  weakly  hydrogen-bonded  materials, 
relaxation  times  conform  to  the  scaling  law 

r  =  ftTvn,  (i) 

where /is  a  function,  V  the  specific  volume,  and  y  a  material 
constant  [4,5].  If  r  is  constant  at  a  transition,  it  follows 
that  the  temperature  Tc  and  volume  Vc  at  the  transition  are 
related  by  the  constancy  of  the  product  variable  Tc  V[ ,  where 
T  is  a  material  constant.  For  a  thermodynamic  transition, 
such  as  for  liquid  crystals,  it  is  not  a  priori  obvious  why 
the  purely  thermodynamic  quantities  Tc  and  Vc  should  be 
related  by  the  exponent  y,  which  scales  the  dynamics  of  the 
system.  However,  the  idea  underlying  density  scaling  is  that  the 
dynamic  properties  of  viscous  liquids  are  governed  primarily 
by  the  repulsive  component  of  the  intermolecular  potential, 
with  changes  in  the  long-range  attractive  part  of  the  potential 
exerting  a  negligible  effect.  To  the  extent  that  this  assumption 
is  valid  for  the  properties  of  interest,  the  potential  can  be 
approximated  by  an  inverse  power  law  (IPL)  or  “soft- sphere” 
intermolecular  potential 

U (r)  oc  r~m .  (2) 

The  IPL  potential  is  scale  invariant:  All  dynamic  properties 
and  excess  thermodynamic  properties  of  the  liquid,  expressed 
in  appropriately  reduced  units,  are  a  function  of  the  scaling 
variable  TVy ,  with  y  =  m/3.  In  particular,  the  relaxation 
times,  viscosities,  and  diffusion  coefficients  (in  reduced  units; 
see  below)  conform  to  the  scaling  law  [Eq.  (1)]  [6,7].  Thus, 
one  can  interpret  the  exponent  y,  derived  from  scaling 
experimental  relaxation  times  or  viscosities,  as  providing  a 
measure  of  the  slope  of  the  effective  intermolecular  potential. 

For  all  state  points  for  which  TVy  is  constant,  not  only  are 
reduced  dynamic  quantities  constant,  but  also  the  structure  of 
the  liquid.  More  generally,  Dyre  and  coworkers  proposed  the 
term  “isomorphs”  for  thermodynamic  paths  along  which  the 
Boltzmann  probabilities  of  scaled  molecular  configurations 
are  invariant  [8,9],  leading  to  the  invariance  of  a  multitude 
of  dynamic  properties.  Liquids  which  have  isomorphs  (the 
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so-called  strongly  correlating  liquids)  also  exhibit  strong 
energy -pres  sure  correlations  [10,1 1]  and,  if  the  proportionality 
constant  for  energy  and  pressure  fluctuations  is  state-point 
independent,  density  scaling.  Liquids  conforming  to  an  IPL  are 
strictly  correlating,  while  those  conforming  to  an  IPL  modified 
with  an  additional  linear  term  are  strongly  correlating  [10,12]. 

For  an  IPL  liquid  the  temperature  and  volume  at  the  melting 
point  are  related  by  the  scaling  law 

TmV£  =  const.,  (3) 

with  r  =  y  =  m/ 3 ;  that  is,  the  melting  curve  represents  an  iso- 
morph  [6].  From  this  follows  the  well-known  Simon  equation 
for  the  pressure  dependence  of  the  melting  temperature  [13] 


where  a ,  To  and  Pq  are  constants,  and  the  parameter  c  = 
(1  +  3 /m)  [14].  Similarly,  reduced  dynamical  quantities  for 
an  IPL  liquid  are  expected  to  be  constant  along  the  melting 
curve  [15]. 

For  a  real  liquid,  the  dynamics  and  the  phase  behavior  are 
expected  to  be  influenced  by  attractive  interactions  and  other 
structure-sensitive  terms  (e.g.,  Coulombic  forces).  However,  a 
large  number  of  materials  (in  particular  nonassociating  molec¬ 
ular  liquids  and  polymers)  exhibit  characteristics  associated 
with  the  IPL  approximation,  including  density  scaling  of  the 
dynamics.  Can  Eq.  (3)  be  generalized  to  the  melting  point  of 
such  liquids?  The  fact  that  liquid  crystal  transitions,  including 
those  from  a  highly  ordered,  crystal-like  mesophase  such  as 
the  smectic  E  to  isotropic  liquid,  conform  to  density  scaling 
suggests  that  this  might  be  the  case.  Dyre  and  coworkers  [8] 
have  predicted  that  along  the  melting  curve  Eq.  (3)  applies  with 
y  =  T,  for  liquids  for  which  the  IPL  is  a  fair  approximation. 
Such  materials  include  nonassociating  liquids  and  liquid 
metals,  but  not,  for  example,  ionic  or  hydrogen-bonded  liquids. 

We  analyze  literature  data  for  a  variety  of  nonassociating 
liquids  to  assess  the  extent  to  which  properties  relating  to 
thermodynamic  transitions,  such  as  density  scaling  of  the 
melting  point  [Eq.  (3)],  can  be  generalized  from  the  idealized 
IPL  to  real  liquids. 

II.  RESULTS 

We  analyze  all  liquids  for  which  the  following  literature 
data  were  available:  (a)  the  dynamics  (viscosity,  self-diffusion 
coefficient,  or  dielectric  relaxation  time)  over  a  substantial 
range  of  temperatures  and  pressures,  with  sufficient  accuracy 
to  permit  superpositioning  in  accord  with  Eq.  (1);  (b)  the 
equation  of  state  (EoS)  over  the  measurement  range  of 
the  dynamic  data  and  in  the  vicinity  of  the  melting  point; 
(c)  the  melting  temperature  as  a  function  of  pressure.  For 
materials  for  which  densities  were  measured  together  with 
viscosities  or  melting  points,  the  former  were  used  rather  than 
an  EoS  from  another  source.  For  liquids  having  more  than  one 
crystalline  form,  we  limited  our  analysis  to  the  lowest-pressure 
crystal  phase.  Crystals  with  a  degree  of  orientational  disorder 
(for  example,  the  low-pressure  crystal  phases  of  methane, 
neopentane,  ethane,  cyclohexane,  carbon  tetrachloride,  and  the 
odd-numbered  n  alkanes  with  n^9)  show  the  same  qualitative 


dependence  of  melting  on  T  and  V  as  orientationally  ordered 
crystals;  thus,  they  were  included  in  this  study. 

A.  Density  scaling  of  the  dynamics 

In  superpositioning  data  according  to  Eq.  (1),  we  use 
reduced  dynamical  quantities,  defined  using  units  of  time 
to  =  vx/3{kT / m)-1/2,  length  Z0  =  u1/3,  and  energy  Eq  =  kT , 
where  mw  and  v  are  molecular  mass  and  volume.  The  reduced 
relaxation  time,  viscosity,  and  diffusion  coefficient  are  then 

r*  =  v~l/3(kT /mw)l/2z, 
r)*  =  i )2/3(mwkT)~l/2ri,  (5) 

D*  =  i rl/3ikT/mwyl/2D. 

The  use  of  reduced  units  is  motivated  by  the  IPL  potential, 
where  density  scaling  applies  to  reduced  quantities  only 
[6,7,15].  The  difference  between  the  scaling  of  reduced 
and  unreduced  quantities  is  negligible  in  the  supercooled 
regime;  however,  most  of  the  liquids  studied  herein  cannot 
be  supercooled,  so  that  the  dynamics  can  only  be  analyzed 
over  a  relatively  narrow  range  above  the  melting  point.  In  such 
cases  the  difference  between  reduced  and  unreduced  quantities 
can  be  substantial  [7];  employing  unreduced  quantities  either 
causes  breakdown  of  the  scaling  or  yields  values  of  the 
exponent  y  which  are  unreasonably  large  or  vary  with  the 
dynamic  variable  (r,  r),  or  D). 

To  date,  density  scaling  has  been  found  to  apply  to  a  large 
number  of  nonassociated  liquids  and  polymers  [4,5]  and  to 
molecular  dynamics  simulations  based  on  IPL  or  Lennard- 
Jones-type  intermolecular  potentials  [9-12,16]  Here  we  show 
that  it  also  holds  for  experimental  viscosities  and  diffusion 
coefficients  of  simple  monoatomic  (argon,  xenon,  krypton)  and 
diatomic  (nitrogen,  oxygen)  liquids,  as  well  as  for  a  variety 
of  other  organic  molecules  for  which  scaling  has  not  been 
previously  reported.  The  obtained  scaling  exponents  are  listed 
in  Table  I  [17-83],  with  representative  plots  of  viscosity  versus 
TVY  displayed  in  Figs.  1  and  2.  For  those  liquids  where  more 
than  one  dynamic  quantity  was  available,  density  scaling  for 
each  gave  exponents  equivalent  within  the  experimental  error. 


B.  Density  scaling  of  the  melting  point 

Figures  3  and  4  are  double  logarithmic  plots  of  the  melting 
temperature  versus  melting  volume  for  a  representative  subset 
of  the  materials  studied.  The  data  can  be  fit  with  straight 
lines.  For  certain  materials,  deviations  from  linearity  were 
observed  either  at  the  highest  or  lowest  pressures,  presumably 
due  to  systematic  errors  in  the  EoS  (which  often  must  be 
extrapolated  to  the  melting  point  using  measurements  at  higher 
temperatures  and/or  lower  pressures).  From  the  linear  fits  we 
define  a  scaling  exponent  for  the  melting  point 


where  T  is  analogous  to  the  well-known  thermodynamic 
potential  parameter  of  liquid  crystal  transitions  [2].  The  range 
of  the  density  and  temperature  data  is  quite  limited,  in  most 
cases  spanning  less  than  one  decade,  as  limited  by  the  extent  of 
the  liquid  phase  in  the  T-V  domain.  For  this  reason,  we  cannot 
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FIG.  1.  (Color  online)  Density  scaling  of  the  reduced  viscosity 
for  selected  aliphatic  and  aromatic  hydrocarbons. 

infer  unambiguously  that  Tm  =  const,  from  the  data  in 
Figs.  3  and  4  alone.  However,  the  parameter  V  is  a  well  defined 
quantity  for  each  liquid,  describing  the  relative  influence  of 
temperature  and  volume  on  the  melting  point.  In  the  case  of 
the  n  alkanes,  the  various  pressure-dependent  melting  points 
were  taken  from  different  studies,  performed  using  different 
techniques,  with  EoS  also  obtained  from  different  sources; 
nevertheless,  the  data  in  Fig.  5  show  a  systematic  increase 
of  the  melting  temperature  and  decrease  of  the  slope  (and 
thus  T)  with  increasing  chain  length.  An  odd-even  effect  is 
also  observed:  odd-numbered  n  alkanes  tend  to  have  higher 


1000/(TVy) 

FIG.  2.  (Color  online)  Density  scaling  of  the  reduced  viscosity 
for  liquid  argon,  krypton,  nitrogen,  carbon  dioxide,  trichloromethane, 
and  carbon  tetrachloride. 


FIG.  3.  (Color  online)  Melting  temperature  vs.  melting  volume 
for  selected  liquids  from  Table  I. 

scaling  exponents  than  the  values  of  even-numbered  ones  in  a 
homologous  series.  A  similar  odd-even  effect  is  well  known  for 
melting  temperatures,  with  even  numbered  n  alkanes  tending 
to  have  the  higher  melting  points  [56]. 

The  same  procedure  can  also  be  applied  to  data  from 
molecular  dynamics  simulations.  Ahmed  et  al.  [84]  reported 
melting  temperatures  and  volumes  for  a  series  of  Lennard- 
Jones  liquids  with  different  repulsive  exponents  ranging  from 
m  =  7  to  12.  The  data  (Fig.  6)  are  well  described  by  the  power 
law  TmVl  —  const,  over  the  entire  range  studied,  with  the 
slope  T  increasing  with  increasing  steepness  of  the  repulsive 
part  of  the  intermolecular  potential  (i.e.,  larger  m). 


C.  Dynamic  versus  thermodynamic  scaling 

In  Fig.  7  we  plot  the  scaling  exponent  for  the  dynamics 
y  versus  the  scaling  exponent  T  for  the  melting  point.  The 


ln(Vm  [cm3/g]) 


FIG.  4.  (Color  online)  Melting  temperature  vs.  melting  volume 
for  the  n- alkanes. 
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FIG.  5.  (Color  online)  Density  scaling  exponents  for  the  melting 
point  and  for  the  dynamics,  for  the  «- alkanes. 


8  - 


6  - 


4  - 


2  - 


1 - > - 1 - * - 1 - 1 - 7| 

/ 

neopentane  / 

•  + 

/ 

/ 

ethane  ' 

x  x 

o,  mo  x 


-  Ar,  Kr,  Xe 


n2,  ch4 


L& 


.  o  o 

*  Dt^ 


'  *■  &<b  H® 
'  B+ 


X  - 


4  6 

Y 


8  10 


diagonal  y  =  T  corresponds  to  viscosity,  relaxation  time, 
and  diffusion  coefficient,  in  reduced  units,  constant  along 
the  melting  curve  (i.e.,  this  curve  is  an  isomorph).  For  eight 
of  the  43  materials  studied,  y  =  T  within  the  experimental 
uncertainty;  these  are  argon,  xenon,  krypton,  nitrogen,  oxygen, 
methane,  ethane,  and  neopentane.  For  the  other  35  materials 
y  >  T,  without  exception. 

The  molecules  of  substances  for  which  y  =  T  are  roughly 
spherical,  rigid,  and  lack  polar  bonds;  this  means  their 
intermolecular  potentials  can  be  approximated  as  spherically 
symmetric.  The  presence  of  polar  bonds  (even  if  the  dipoles 
cancel  out  so  there  is  no  net  dipole  moment)  breaks  the 
symmetry  of  the  potential.  Thus,  although  y  =  T  for  methane, 
y  >  T  for  carbon  tetrachloride  and  tetrafluoride,  which  have 


FIG.  7.  (Color  online)  Scaling  exponent  for  the  melting  point 
vs.  for  the  dynamics,  for  the  liquids  in  Table  I:  nonpolar  spher¬ 
ical  molecules  (circles),  n  alkanes  (squares),  branched  and  cyclic 
alkanes  (crosses),  aromatic  hydrocarbons  (hexagons),  halomethanes 
(diamonds),  various  other  molecules:  polar  (triangles)  and  nonpolar 
(crosses). 

the  same  molecular  configuration  as  methane  but  polar  carbon- 
halogen  bonds.  The  same  observation  applies  to  isopentane 
(y  =  T)  compared  to  tetramethylsilane  (y  >  T),  where  the 
C-C  bonds  of  the  former  are  replaced  by  (weakly)  polar  C-S 
bonds. 

Similarly,  departure  from  a  quasispherical  shape  results  in 
unequal  melting  and  dynamic  scaling  exponents.  For  example, 


InV 

m 

FIG.  6.  (Color  online)  Melting  temperature  vs.  melting  volume 
for  the  m- 6  Lennard- Jones  system  (data  from  Ref.  [84])  The  slope 
(scaling  exponent  for  the  melting  point)  increases  with  increasing 
steepness  m  of  the  repulsive  part  of  the  potential. 


FIG.  8.  (Color  online)  Difference  between  scaling  exponents  for 
the  dynamics  (y)  and  the  melting  point  (r)  as  a  function  of  the 
acentric  factor  go.  Symbols  are  as  in  Fig.  7. 
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TABLE  I.  Density  scaling  exponent  for  the  melting  point  (T)  and  for  dynamics  (y),  acentric  factor  (from  Ref.  [17]),  and  literature 
sources  for  dynamics,  melting  points  and  equation  of  state  for  the  liquids  studied.  Temperature  ranges  and  maximum  pressures  for  the 
dynamical  data,  as  well  as  the  range  of  melting  points  used  in  the  scaling  procedures,  are  also  given. 


No. 

Substance 

r 

Y 

co 

T  [K] 

/’max  [MPa] 

Tm  [K] 

t]  or  D  or  z(T,P) 

TJP) 

V(T,P) 

1 

argon 

4.0 

4.0 

0.000 

90-500 

400 

84-348 

[18]  ft) 

[19,20] 

[21] 

2 

krypton 

4.0 

4.0 

0.000 

298-348 

200 

116-157 

[21]  ft) 

[22] 

[23] 

3 

xenon 

4.1 

4.0 

0.000 

248-343 

135 

161-289 

[24]  (D) 

[25] 

[18] 

4 

nitrogen 

3.8 

4.0 

0.040 

285-473 

2200 

63-191 

[18]  ft) 

[26] 

[26] 

5 

oxygen 

4.8 

5.0 

0.022 

54-200 

82 

55-63 

[18]  ft) 

[27] 

[27] 

6 

carbon  dioxide 

1.9 

4.5 

0.228 

220-300 

453 

216-265 

[28]  ft) 

[29] 

[29] 

7 

dichloromethane 

4.2 

5.0 

0.192 

186-306 

200 

177-200 

[30]  (D) 

[31] 

[32] 

8 

carbon  tetrachloride 

4.9 

7.0 

0.193 

283-328 

148 

289-400 

[33]  (D) 

[34] 

[32] 

9 

carbon  tetrafluoride 

4.9 

6.0 

0.186 

140-432 

200 

89-224 

[35]  (D) 

[36] 

[18] 

10 

acetonitrile 

2.6 

3.4 

0.338 

298-373 

500 

229-300 

[37]  ft) 

[38] 

[39] 

11 

pyridine 

2.3 

4.2 

0.239 

303-423 

500 

334-346 

[40]  (D,rj) 

[41] 

[39] 

12 

salol 

3.7 

5.3 

0.632 

265-380 

700 

314-383 

[42]  (r) 

[38] 

[43] 

13 

octamethyltetrasiloxane 

8.4 

9.5 

0.589 

293-493 

200 

289-363 

[44]  (D) 

[44] 

[39] 

14 

hexafluorobenzene 

3.6 

5.5 

0.395 

288-493 

200 

277-343 

[45](D) 

[46] 

[32] 

15 

tetramethylsilane 

4.5 

5.5 

0.175 

298-373 

455 

165-224 

[47]  ft) 

[48] 

[32] 

16 

methane 

4.2 

4.0 

0.011 

95-180 

250 

91-255 

[18]  ft) 

[49] 

[49] 

17 

ethane 

5.9 

5.0 

0.099 

136-454 

200 

90-185 

[50]  (D) 

[51] 

[51] 

18 

propane 

4.1 

7.0 

0.152 

112-453 

200 

86-125 

[50]  (D) 

[52] 

[52] 

19 

?z-butane 

3.2 

5.8 

0.199 

135-400 

69 

113-125 

[18]  ft) 

[53] 

[53] 

20 

n- pentane 

3.4 

7.5 

0.249 

303-450 

250 

143-184 

[54]  ft),  [55]  ft) 

[56] 

[57,58] 

21 

/7-hexane 

2.9 

7.0 

0.305 

298-450 

500 

177-225 

[55]  ft),  [59]  ft) 

[56] 

[57,58] 

22 

n- heptane 

2.6 

6.3 

0.351 

303-323 

69 

182-225 

[60]  ft) 

[56] 

[57,58] 

23 

/z-octane 

2.7 

6.2 

0.396 

283-473 

375 

216-269 

[61]  ft),  [62]  ft) 

[56] 

[57,58] 

24 

/z-nonane 

2.3 

5.9 

0.438 

303-323 

69 

220-268 

[60]  ft) 

[56] 

[57,58] 

25 

/z-decane 

2.1 

5.9 

0.484 

298-373 

250 

243-294 

[55]  ft),  [62]  ft) 

[56] 

[57,58] 

26 

/z-undecane 

1.7 

5.3 

0.536 

303-323 

69 

248-295 

[60]  ft) 

[63] 

[64] 

27 

n- dodecane 

2.0 

5.2 

0.573 

298-473 

200 

264-322 

[65]  ft) 

[63] 

[64] 

28 

/z-hexadecane 

1.9 

4.7 

0.747 

298-373 

425 

291-348 

[59]  ft) 

[66] 

[66] 

29 

/z-octadecane 

1.8 

4.9 

0.795 

323-473 

90 

301-353 

[65]  ft) 

[63] 

[67] 

30 

neopentane 

8.2 

8.4 

0.196 

267-450 

120 

257-390 

[68]  (D) 

[68] 

[69] 

31 

isopentane 

5.1 

9.0 

0.228 

298-328 

200 

113-151 

[70]  (D) 

[56] 

[69] 

32 

cyclopentane 

6.4 

9.0 

0.194 

298-328 

200 

178-251 

[70]  (D) 

[56] 

[69] 

33 

cyclohexane 

6.3 

7.5 

0.212 

313-383 

214 

279-417 

[71]  (D.r,) 

[56] 

[69] 

34 

methylcyclohexane 

5.8 

7.3 

0.235 

298-363 

500 

147-196 

[72]  (D),  [73]  ft) 

[56] 

[69] 

35 

benzene 

2.7 

5.0 

0.211 

303-433 

400 

306-437 

[74]  ft) 

[75] 

[76] 

36 

toluene 

4.5 

5.7 

0.264 

298-363 

149 

177-900 

[61]  ft),  [73]  ft) 

[77] 

[76] 

37 

oxylene 

3.2 

7.4 

0.313 

298-363 

100 

248-298 

[73]  ft) 

[38] 

[76] 

38 

m-xylene 

3.4 

7.5 

0.326 

298-473 

200 

228-283 

[62]  ft),  [73]  ft) 

[38] 

[76] 

39 

p-xylene 

3.1 

7.2 

0.326 

313-363 

100 

287-337 

[73]  ft) 

[38] 

[76] 

40 

mesitylene 

5.9 

6.5 

0.398 

298-313 

280 

221-306 

[78]  (D) 

[78] 

[76] 

41 

1  -methylnaphthalene 

2.0 

6.3 

0.348 

298-473 

200 

242-314 

[62]  ft) 

[75] 

[76] 

42 

o-terphenyl 

3.5 

4.0 

0.467 

256-308 

125 

329-477 

[79]  (r),  [80]  (r) 

[81] 

[82] 

43 

/z-dodecylbenzene 

2.5 

4.8 

0.786 

298-363 

100 

269-292 

[73]  ft) 

[83] 

[73] 

in  the  case  of  the  /i-alkanes,  the  roughly  spherical  methane  and 
ethane  have  y  =  T,  but  for  longer  chains  (having  an  elongated 
molecular  shape  and  intramolecular  degrees  of  freedom) 

y  >  r. 

A  parameter  widely  used  to  quantify  the  deviation  of  a 
molecule  from  sphericity  is  the  acentric  factor  co,  introduced 
by  Pitzer  [85].  This  factor,  calculated  from  the  critical 
properties  of  the  liquid,  is  defined  as  co  =  —  log(Pr)  —  1  at 
Tr  =0.7,  where  P,  and  Tr  are  the  reduced  vapor  pressure 
and  temperature  relative  to  the  critical  point.  A  value  of  co  =  0 
corresponds  to  a  perfectly  spherical  molecule.  Acentric  factors 


for  the  liquids  studied  are  listed  in  Table  I;  they  vary  from  0 
to  0.8.  In  Fig.  8,  the  difference  y  —  V  is  plotted  as  a  function 
of  co.  Only  those  liquids  for  which  y  =  F  have  co  <  0.1,  with 
the  exception  of  neopentane,  for  which  co  =  0.196.  For  the 
remaining  liquids  co  >  0.1  and  y  >  T,  although  there  is  no 
correlation  between  the  magnitude  the  acentric  factor  and  the 
difference  between  the  two  scaling  exponents. 

We  also  note  there  is  no  correlation  of  the  difference 
between  y  and  V  with  properties  such  as  the  molecular  dipole 
moment,  or  the  excess  enthalpy  or  entropy  of  melting.  The 
crystal  structure  also  appears  to  be  irrelevant:  Of  the  materials 


031504-5 


D.  FRAGIADAKIS  AND  C.  M.  ROLAND 


PHYSICAL  REVIEW  E  83,  031504  (2011) 


conforming  to  y  =  T,  argon,  krypton,  xenon,  methane,  and 
neopentane  form  face-centered  cubic  crystals,  while  the  low 
pressure  crystal  phase  of  nitrogen  is  hexagonal,  oxygen  is 
monoclinic,  and  ethane  is  body-centered  cubic.  Moreover, 
the  low-pressure  crystal  phases  of  methane,  neopentane, 
and  ethane  are  plastic  crystals,  possessing  some  degree  of 
orientational  disorder.  However,  plastic  crystallinity  does  not 
give  rise  to  equivalent  values  for  y  and  T :  Carbon  tetrafluoride 
and  carbon  tetrachloride,  tetramethylsilane,  and  certain  odd- 
numbered  alkanes  have  orientationally  disordered  crystals,  but 
we  find  y  >  T . 

As  stated,  y  is  always  greater  than  or  equal  to  T;  in  no 
case  do  we  observe  T  >  y .  This  implies  that  deviation  from 
spherical  symmetry  of  the  intermolecular  potential  causes  the 
dynamics  to  be  more  sensitive  to  volume  than  the  melting 
point.  This  means  that  the  effective  intermolecular  potential 
for  dynamical  properties  is  always  steeper  than  that  relevant  to 
crystal  melting.  Whereas  the  dynamics  reflect  an  orientational 
average  of  the  intermolecular  potential,  the  melting  point 
is  affected  predominantly  by  the  potential  for  molecules  at 
specific  angles,  corresponding  to  their  relative  orientation  in 
the  crystal  lattice.  Moreover,  as  volume  is  reduced,  molecular 
jamming  increases  the  viscosity,  with  these  strong  steric 
interactions  exacerbated  for  molecules  having  nonspherical 
shapes  [86].  However,  such  effects  do  not  obviously  influence 
the  relative  stability  of  the  crystal  phase. 


III.  DISCUSSION 

A.  Scaling  exponents  and  molecular  structure 

Collected  in  Table  I  are  density  scaling  exponents  for  43 
liquids  having  relatively  simple  molecular  structures.  From 
these  data  we  can  draw  inferences  concerning  the  relationship 
between  y ,  which  is  governed  by  the  steepness  of  the  effective 
intermolecular  potential,  and  structural  properties.  Polar  liq¬ 
uids  tend  to  have  low  values  of  y ,  in  agreement  with  both  prior 
experimental  results  [87]  and  molecular  dynamics  simulations 
in  which  adding  a  permanent  dipole  to  an  asymmetric  two- site 
Lennard-Jones  molecule  caused  a  similar  decrease  in  y  [9]. 

Increasing  molecular  flexibility  (i.e.,  more  internal  degrees 
of  freedom)  corresponds  to  lower  values  of  y.  For  example, 
the  ^-alkanes  have  smaller  scaling  exponents  than  the  more 
rigid  branched  and  cyclic  alkanes.  This  reflects  the  softening 
of  the  interatomic  potential  due  to  the  pressure-insensitive 
intramolecular  bonds  [87,88].  For  the  homologous  series  of 
/i-alkanes  having  moderate  (pentyl)  or  long  chain  lengths,  both 
the  dynamic  and  melting-point  scaling  exponents  continuously 
decrease  with  increasing  chain  length  (Fig.  5),  with  the 
difference  between  the  two  remaining  approximately  constant. 

There  are  some  cases  for  which  a  systematic  change  in 
chemical  structure  does  not  lead  to  corresponding  changes 
in  scaling  exponents.  For  example,  the  presence  of  pendant 
methyl  groups  on  benzene  does  not  lead  to  systematic  changes 
in  y,  as  seen  by  comparing  data  for  toluene,  xylenes,  and 
mesitylene.  In  this  respect  the  scaling  behavior  is  similar 
to  other  properties  (viscosity,  glass  transition  temperature, 
nonexponentiality  of  relaxation,  fragility,  melting  point,  and 
thermodynamic  properties  such  as  melting  and  other  transi¬ 
tions),  which  likewise  do  not  correlate  in  any  obvious  way 


with  structure  (that  is,  these  fundamental  properties  cannot 
be  determined  a  priori  from  knowledge  of  the  chemical 
structure).  Subtle  effects  related  to  packing  and  symmetries 
exert  a  nontrivial  influence. 


B.  Relationship  to  thermodynamic  Griineisen  parameter 

According  to  the  Lindemann  criterion,  melting  transpires 
when  the  mean  vibrational  amplitude  reaches  a  critical  fraction 
of  the  interatomic  distance  [14].  Gilvarry  [89]  derived  a  form 
of  the  Lindemann  equation  (the  so-called  Lindemann- Gilvarry 
relation,  widely  used  in  geophysics  [90])  that  relates  the  change 
in  volume  along  the  melting  curve  to  the  thermodynamic 
Griineisen  parameter 


T 


d\nTm 
d  In  Vm 


=  2 (yG  -  |). 


(7) 


The  dynamic  scaling  exponent  can  also  be  related  to  the 
Griineisen  parameter,  based  on  the  IPL  approximation  [91] 


Y  =  2 (Vg  ~  |),  (8) 

or  from  the  assumption  that  the  dynamics  are  governed  by  the 
excess  entropy  [92,93] 


Y  = 


Cy 

ACi 


~Yg- 


(9) 


where  the  prefactor  is  the  inverse  of  the  relative  difference  in 
isochoric  heat  capacity  for  the  liquid  and  crystal.  A  related 
formula  is  [12] 


(10) 


Equations  (7)— (10)  indicate  a  link  between  the  dynamics 
and  the  melting  point  that  may  extend  more  broadly  to  liquids 
for  which  an  IPL  potential  is  not  applicable.  However,  since 
empirically  y  ~  T  is  observed  only  for  a  subset  of  the  liquids 
in  Table  I,  modification  of  the  Lindemann  criterion  or  the 
arguments  leading  to  Eqs.  (8)-(10)  is  necessary  for  consistency 
with  experimental  results  for  nonspherical  liquids  and  those 
possessing  internal  degrees  of  freedom. 

Molecular  dynamics  simulations  have  been  used  to  study 
the  effect  of  molecular  shape  on  the  scaling  exponent.  Using 
a  repulsive  exponent  m  =  12,  a  density  scaling  exponent  of 
y  =  5.0  was  obtained  for  a  binary  mixture  of  Lennard-Jones 
spheres  [16].  For  an  asymmetric  dumbbell  molecule  consisting 
of  two  different  sized  Lennard-Jones  spheres  (also  with  m  = 
12),  the  scaling  exponent  increases  to  6.1.  And  for  the  Lewis- 
Wahnstrom  ortho-terphenyl  (OTP)  model,  consisting  of  three 
rigidly  connected  spheres,  y  =  7.9  [9].  Melting  point  data 
were  not  reported,  but  it  would  be  of  interest  to  determine 
T  for  these  asymmetric  molecules.  The  expectation  from  the 
behavior  of  real  materials  is  that  T  <  y . 


C.  Viscosity  at  the  melting  point 

The  dependence  of  liquid  dynamics  on  pressure  and  temper¬ 
ature  is  important  in  the  field  of  geophysics.  Characterization 
of  the  relevant  materials  (iron  and  its  alloys,  silicates)  is 
difficult  because  melting  occurs  at  extreme  temperatures  and 
pressures  of  tens  of  GPa  or  more.  Thus,  extrapolations  are 
required  of  properties  measured  for  experimentally  attainable 
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conditions.  One  approach  is  to  assume  that  the  viscosity 
is  constant  at  the  melting  point,  independent  of  pressure. 
However,  experimental  data  on  geological  materials  indicate 
inferentially  an  increase  of  r\  along  the  melting  line  with 
increasing  P  [94,95].  This  implies  that  at  sufficiently  high  pres¬ 
sures,  geophysical  materials  might  vitrify  at  the  melting  point, 
suggesting  that  parts  of  the  earth’s  core  [94]  or  lower  mantle 
[95]  are  in  the  glassy  state.  Of  course,  such  conclusions  are 
speculative  given  the  absence  of  direct  measurements  [96,97]. 

The  extension  of  the  results  herein  for  molecular  liquids  to 
materials  of  geological  interest  could  clarify  this  issue.  Note 
that  liquid  metals  are  predicted  to  be  strongly  correlating  [8], 
and  accordingly  density  scaling  should  apply.  Since  metals 
are  monoatomic,  we  anticipate  that  y  =  T.  If  this  is  the  case, 
dynamical  quantities  at  the  melting  point  are  constant  when 
expressed  in  reduced  quantities  so  that  the  viscosity  of  metals 
at  the  melting  point  can  be  expressed  as 

t){Tm)  cx  V~2/3T^2.  (11) 

This  equation  (with  a  universal  prefactor)  was  proposed 
empirically  by  Andrade  [98,99]  (in  units  of  poise): 

ri(Tm)  =  5.1  x  10-4(A^/2,  (12) 

where  A  is  the  atomic  weight  and  Va  the  molar  volume  at  the 
melting  point.  Equation  (12)  has  been  used  to  extrapolate  the 


viscosity  at  ambient  pressure  of  various  metals  to  the  melting 
point.  In  liquids  for  which  y  >  T,  the  viscosity  at  the  melting 
temperature  would  increase  more  strongly  with  pressure. 

IV.  SUMMARY 

Two  scaling  relations,  for  the  dynamics  and  for  the 
melting  behavior,  yield  respective  exponents  y  and  T.  For 
spherical  molecules  lacking  internal  dipole  moments,  we 
find  these  exponents  to  be  equivalent;  however,  more  gen¬ 
erally  y  >  T.  The  former  condition  prevails  exactly  for  an 
IPL  intermolecular  potential,  reflecting  the  appropriateness 
of  the  IPL  approximation  for  many  materials.  From  the 
scaling  behavior  two  empirical  relations,  the  Lindemann- 
Gilvarry  and  the  Andrade  equations,  can  be  derived.  In 
addition  to  providing  insights  regarding  the  forces  between 
molecules,  the  scaling  exponents  can  be  used  to  estimate 
melting  points  at  experimentally  inaccessible  temperatures  and 
pressures. 

ACKNOWLEDGMENTS 

This  work  was  supported  by  the  Office  of  Naval  Research. 
We  thank  R.  Casalini  for  useful  discussions.  D.F.  acknowl¬ 
edges  the  National  Research  Council  for  a  postdoctoral 
fellowship. 


[1]  C.M.  Roland,  Soft  Matter  4,  2316  (2008). 

[2]  S.  Urban  and  C.  M.  Roland,  J.  Non-Cryst.  Solids  357,  740 
(2010). 

[3]  C.  M.  Roland,  R.  B.  Bogoslovov,  R.  Casalini,  A.  R.  Ellis,  S.  Bair, 
S.  J.  Rzoska,  K.  Czuprynski,  and  S.  Urban,  J.  Chem.  Phys.  128, 
224506  (2008). 

[4]  C.  M.  Roland,  S.  Hensel-Bielowka,  M.  Paluch,  and  R.  Casalini, 
Rep.  Prog.  Phys.  68,  1405  (2005). 

[5]  C.  M.  Roland,  Macromolecules  43,  7875  (2010). 

[6]  W.  G.  Hoover  and  M.  Ross,  Contemp.  Phys.  12,  339 
(1971). 

[7]  D.  Fragiadakis  and  C.  M.  Roland,  J.  Chem.  Phys.  134,  044504 
(2011). 

[8]  N.  Gnan,  T.  B.  Schrpder,  U.  R.  Pedersen,  N.  P.  Bailey,  and  J.  C. 
Dyre,  J.  Chem.  Phys.  131,  234504  (2009). 

[9]  T.  B.  Schrpder,  U.  R.  Pedersen,  N.  P.  Bailey,  S.  Toxvaerd,  and 
J.  C.  Dyre,  Phys.  Rev.  E  80,  041502  (2009). 

[10]  N.  P.  Bailey,  U.  R.  Pedersen,  N.  Gnan,  T.  B. 
Schrpder,  and  J.  C.  Dyre,  J.  Chem.  Phys.  129,  184507 
(2008);  129,  184508  (2008). 

[11]  D.  Coslovich  and  C.  M.  Roland,  J.  Chem.  Phys.  130,  014508 
(2009). 

[12]  T.  B.  Schrpder,  N.  P.  Bailey,  U.  R.  Pedersen,  N.  Gnan,  and  J.  C. 
Dyre,  J.  Chem.  Phys.  131,  234503  (2009). 

[13]  F.  E.  Simon  and  G.  Glatzel,  Z.  Anorg.  Allg.  Chem.  178,  300 
(1929). 

[14]  N.  H.  March  and  M.  P.  Tosi,  Introduction  to  Liquid  State  Physics 
(World  Scientific,  Singapore,  2002). 

[15]  Y.  Hiwatari,  H.  Matsuda,  T.  Ogawa,  N.  Ogita,  and  A.  Ueda, 
Prog.  Theor.  Phys.  52,  1105  (1974). 


[16]  D.  Coslovich  and  C.  M.  Roland,  J.  Phys.  Chem.  B  112,  1329 
(2008). 

[17]  Chemical  Properties  Handbook ,  edited  by  C.  L.  Yaws 
(McGraw-Hill,  New  York,  1999). 

[18]  NIST  Chemistry  WebBook,  NIST  Standard  Reference  Database 
Number  69,  edited  by  P.  J.  Linstrom  and  W.  G.  Mallard  (NIST, 
Gaithersburg,  MD,  2005). 

[19]  Ch.  Tegeler,  R.  Span,  and  W.  Wagner,  J.  Phys.  Chem.  Ref.  Data 
28,  779  (1999). 

[20]  S.  M.  Stishov,  Sov.  Phys.  Usp.  17,  625  (1975). 

[21]  N.  J.  Trappeniers,  A.  Botzen,  J.  Van  Oosten,  and  H.  R.  van  der 
Berg,  Physica  31,  945  (1965). 

[22]  A.  Michels  and  C.  Prins,  Physica  28,  101  (1962). 

[23]  W.  B.  Streett  and  L.  A.  K.  Staveley,  J.  Chem.  Phys.  55,  2495 
(1971). 

[24]  P.  W.  E.  Peereboom,  H.  Luigjes,  and  K.  O.  Prins,  Physica  A 156, 
260  (1989). 

[25]  W.  B.  Streett,  L.  S.  Sagan,  and  L.  A.  K.  Steveley,  J.  Chem. 
Thermodyn.  5,  633  (1973). 

[26]  R.  Span,  E.  W.  Lemmon,  R.  T.  Jacobsen,  W.  Wagner,  and 
A.  Yokozeki,  J.  Phys.  Chem.  Ref.  Data  29,  1361  (2000). 

[27]  R.  B.  Stewart,  R.  T.  Jacobsen,  and  W.  Wagner,  J.  Phys.  Chem. 
Ref.  Data  20,917  (1991). 

[28]  P.  S.  van  der  Gulik,  Physica  A  238,  81  (1997). 

[29]  R.  Span  and  W.  Wagner,  J.  Phys.  Chem.  Ref.  Data  25,  1509 
(1996). 

[30]  F.  X.  Prielmeier  and  H.-D.  Liidemann,  Mol.  Phys.  58,  593 
(1986). 

[31]  E.  Yu.  Tonkov,  High  Pressure  Phase  Transformations,  A 
Handbook ,  Vol.  1  (Gordon  &  Breach,  Philadelphia,  1992). 


031504-7 


D.  FRAGIADAKIS  AND  C.  M.  ROLAND 


PHYSICAL  REVIEW  E  83,  031504  (2011) 


[32]  I.  Cibulka,  T.  Takagi,  and  K.  Rzicka,  J.  Chem.  Eng.  Data  46,  2 

(2001). 

[33]  M.  A.  McCool  and  L.  A.  Woolf,  J.  Chem.  Soc.,  Faraday  Trans. 
1  68,  1971  (1972). 

[34]  R.  J.  Zhu,  W.  Xu,  Y.  L.  Tian,  and  J.  S.  Hao,  Chin.  Phys.  B  17, 
1088  (2008). 

[35]  M.  Has  and  H.-D.  Liidemann,  J.  Mol.  Liq.  46,  7  (1990). 

[36]  D.  van  der  Putten,  K.  O.  Prins,  P.  J.  Kortbeek,  and  N.  J. 
Trappeniers,  J.  Phys.  C:  Solid  State  Phys.  20,  3161  (1987). 

[37]  J.  H.  Dymond,  M.  A.  Awan,  N.  F.  Glen,  and  J.  D.  Isdale,  Int.  J. 
Thermophysics  12,  433  (1991). 

[38]  S.  E.  Babb,  Rev.  Mod.  Phys.  35,  400  (1963). 

[39]  I.  Cibulka  and  T.  Takagi,  J.  Chem.  Eng.  Data  47,  1037  (2002). 

[40]  M.  Fury,  G.  Munie,  and  J.  Jonas,  J.  Chem.  Phys.  70, 1260  (1979). 

[41]  T.  Ebina,  C.  Yokohama,  and  S.  Takahashi,  J.  Jap.  Petroleum  Inst. 
36,315  (1993). 

[42]  R.  Casalini,  C.  M.  Roland,  and  M.  Paluch,  J.  Phys.  Chem.  A 
107,  2369  (2003). 

[43]  L.  Comez,  S.  Corezzi,  D.  Fioretto,  H.  Kriegs,  A.  Best,  and 
W.  Steffen,  Phys.  Rev.  E  70,  011504  (2004). 

[44]  A.  Greiner- Schmid,  M.  Has,  and  H.-D.  Liidemann,  Z. 
Naturforsch.  Sect.  A-J.  Phys.  Sci.  45a,  1281  (1990). 

[45]  D.  L.  Hogenboom,  K.  Krynicki,  and  D.  W.  Sawyer,  Mol.  Phys. 
40,  823  (1980). 

[46]  M.  Z.  Faizullin,  Russian  J.  Phys.  Chem.  59,  245  (1985). 

[47]  H.  J.  Parkhurst  Jr.  and  J.  Jonas,  J.  Chem.  Phys.  63,  2705  (1975). 

[48]  T.  Atake  and  H.  Chihara,  Chem.  Phys.  Lett.  56,  330  (1978). 

[49]  U.  Setzmann  and  W.  Wagner,  J.  Phys.  Chem.  Ref.  Data  20, 1061 
(1991). 

[50]  A.  Greiner-Schmid,  S.  Wappmann,  M.  Has,  and  H.-D. 
Liidemann,  J.  Chem.  Phys.  94,  5643  (1991). 

[51]  D.  Biicker  and  W.  Wagner,  J.  Phys.  Chem.  Ref.  Data  35,  205 
(2006). 

[52]  E.  W.  Lemmon,  M.  O.  McLinden,  and  W.  Wagner,  J.  Chem. 
Eng.  Data  54,3141  (2009). 

[53]  D.  Biicker  and  W.  Wagner,  J.  Phys.  Chem.  Ref.  Data  35,  929 
(2006). 

[54]  E.  Kiran  and  Y.  L.  Sen,  Int.  J.  Thermophys.  13,  411  (1992). 

[55]  C.  M.  B.  P.  Oliveira  and  W.  A.  Wakeham,  Int.  J.  Thermophys. 
13,  773  (1992). 

[56]  A.  Wiirflinger,  Berichte  der  Bunsengesellschaft  fur  Physikalis- 
che  Chemie  79,  1195  (1975). 

[57]  I.  Cibulka,  Fluid  Phase  Equilib.  89,  1  (1993). 

[58]  I.  Cibulka  and  L.  Hnedkovsky,  J.  Chem.  Eng.  Data  41,  657 
(1996). 

[59]  J.  H.  Dymond,  K.  J.  Young,  and  J.  D.  Isdale,  Int.  J.  Thermophys. 

I,  345  (1980). 

[60]  M.  J.  Assael,  M.  Papadaki,  and  W.  A.  Wakeham,  Int.  J. 
Thermophys.  12,  449  (1991). 

[61]  K.  R.  Harris,  R.  Malhotra,  and  L.  A.  Woolf,  J.  Chem.  Eng.  Data 
42,  1254  (1997). 

[62]  D.  R.  Caudwell,  J.  P.  M.  Trusler,  V.  Vesovic,  and  W.  A.  Wakeham, 

J.  Chem.  Eng.  Data  54,  359  (2009). 

[63]  A.  Wiirflinger  and  G.  M.  Schneider,  Berichte  der  Bunsenge¬ 
sellschaft  fur  Physikalische  Chemie  77,  121  (1973). 

[64]  R.  Landau  and  A.  Wiirflinger,  Berichte  der  Bunsengesellschaft 
fur  Physikalische  Chemie  84,  895  (1980). 


[65]  D.  R.  Caudwell,  J.  P.  M.  Trusler,  V.  Vesovic,  and  W.  A.  Wakeham, 
Z.  Naturfors.  Sect.  A-J.  Phys.  Sci.  25,  1339  (2004). 

[66]  A.  Wiirflinger  and  M.  Sandmann,  Z.  Naturfors.  Sect.  A-J.  Phys. 
Sci.  55,  533  (2000). 

[67]  A.  Wiirflinger,  D.  Mondieig,  F.  Rajabalee,  and  M.  A.  Cuervas- 
Diarte,  Z.  Naturfors.  Sect.  A-J.  Phys.  Sci.  56,  626  (2001). 

[68]  M.  Woznyj,  F.  X.  Prielmeier,  and  H.-D.  Liidemann,  Z. 
Naturforsch.  A  39,  800  (1984). 

[69]  I.  Cibulka  and  T.  Takagi,  J.  Chem.  Eng.  Data  44,  1105  (1999). 

[70]  A.  Enninghorst,  F.  D.  Wayne,  and  M.  D.  Zeidler,  Mol.  Phys.  88, 
437  (1996). 

[71]  J.  Jonas,  D.  Hasha,  and  S.  G.  Huang,  J.  Chem.  Phys.  84,  109 
(1980). 

[72]  J.  Jonas,  D.  Hasha,  and  S.  G.  Huang,  J.  Chem.  Phys.  71,  3996 
(1979). 

[73]  A.  Et-Tahir,  C.  Boned,  B.  Lagourette,  and  P.  Xans,  Int.  J. 
Thermophys.  16,  1309  (1995). 

[74]  H.  J.  Parkhurst  Jr.  and  J.  Jonas,  J.  Chem.  Phys.  63,  2705  (1975). 

[75]  C.  Yokoyama,  T.  Ebina,  and  S.  Takahashi,  Fluid  Phase  Equilib. 
84,  207  (1993). 

[76]  I.  Cibulka  and  T.  Takagi,  J.  Chem.  Eng.  Data  44,  411  (1999). 

[77]  J.  Osugi,  K.  Shimizu,  K.  Yasunami,  M.  Moritoki,  and  A. 
Onodera,  Rev.  Phys.  Chem.  Japan  38,  90  (1968). 

[78]  A.  J.  Easteal  and  L.A.  Woolf,  Int.  J.  Thermophys.  8,  71  (1987). 

[79]  M.  Naoki,  H.  Endou,  and  K.  Matsumoto,  J.  Chem.  Phys.  91, 
4169  (1987). 

[80]  G.  Fytas,  Th.  Dorfmiiller,  and  C.  H.  Wang,  J.  Chem.  Phys.  87, 
5041  (1983). 

[81]  V.  J.  Fratello  and  V.  E.  Bean,  Int.  J.  Thermophys.  4,  253 
(1983). 

[82]  M.  Naoki  and  S.  Koeda,  J.  Phys.  Chem.  93,  948  (1989). 

[83]  M.  Milhet,  J.  Pauly,  J.  A.  P.  Coutinho,  and  J.-L.  Daridon, 
J.  Chem.  Eng.  Data  53,  233  (2008). 

[84]  A.  Ahmed  and  R.  J.  Sadus,  J.  Chem.  Phys.  131,  174504  (2009). 

[85]  J.  H.  Robson,  J.  Am.  Chem.  Soc.  77,  107  (1955). 

[86]  K.  L.  Ngai  and  C.  M.  Roland,  Macromolecules  26,  6824  (1993). 

[87]  C.  M.  Roland,  S.  Bair,  and  R.  Casalini,  J.  Chem.  Phys.  125, 
124508  (2006). 

[88]  R.  Casalini,  C.  M.  Roland,  and  S.  Capaccioli,  J.  Chem.  Phys. 
126,  184903  (2007). 

[89]  J.  J.  Gilvarry,  Phys.  Rev.  102,  308  (1956). 

[90]  J.  P.  Poirier,  Introduction  to  the  Physics  of  the  Earth’s  Interior 
(Cambridge  University  Press,  Cambridge,  England,  2000). 

[91]  C.  M.  Roland,  J.  L.  Feldman,  and  R.  Casalini,  J.  Non-Cryst. 
Solids  352,  4895  (2006). 

[92]  R.  Casalini,  U.  Mohanty,  and  C.  M.  Roland,  J.  Chem.  Phys.  125, 
014505  (2006). 

[93]  C.  M.  Roland  and  R.  Casalini,  J.  Phys.  Condens.  Matter  19, 
205118  (2007). 

[94]  V.  V.  Brazhkin  and  A.  G.  Lyapin,  Phys.  Usp.  43,  493  (2000); 
M.  Z.  Faizullin  and  V.  P.  Skripov,  High  Temp.  40,  364  (2002). 

[95]  I.  Avramov,  J.  Phys.  Condens.  Matter  20,  244106  (2008). 

[96]  V.  V.  Brazhkin,  J.  Phys.  Condens.  Matter  20,  244102  (2008). 

[97]  V.  N.  Mineev  and  A.  I.  Funtikov,  Phys.  Usp.  47,  671 
(2004). 

[98]  E.  N.  da  C.  Andrade,  Nature  128,  835  (1931). 

[99]  E.  N.  da  C.  Andrade,  Philos.  Mag.  Series  7,  112,  497  (1934). 


031504-8 


